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The nature of this project was unfortunately, and unrecognisabr ^d?/ 
foreshadowed by events bounding ay participation in the first Landsat-D 
Investigators Workshop: aircraft deviation around a tornado and unscheduled 

«rv refueling resulted in landing at Dallas Ft. Worth Airport 30 minutes after 

Branlff declared bankrupty, causing, amongest other complications, missing the 
4 0 first half day of presentations; materials gathered at the workshop and mailed 
at the primary U.S. Postal Office (Union Station) never arrived at my home 
^ office. With the best intentions of many persons, the caprice and perversity 
^ of this inantlmate behemoth reluctantly yields. 
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Our work was from the onset intended to be a detailed examination of the 
geometric and radiometric properties of the Landsat 4 imaging data , in 
contrast to an application analyses of the fully processed Images. 
Consequently, our request for data and information have been more extensive 
and less well specified than probably anticipated. We initially specified 
some of the required data in terms of their radiometric content, rather than 
by specific location. This appears not to have been effective; and we will 
respecify these data requirements in terms of geographic and temporal 
constraints. 

The lack of information from the project in terms of data availablllity 
and its plans for providing data to us has been frustrating. We have thus far 
received data for one of the scenes specified in our data requirements 
(Washington, D.C.), and data tapes for two other scenes used primarily for 
checkiug our ability to handle the Thematic Mapper formats. It has required 
repeated personal requests to obtain a list of the Thematl: Mapper scenes 
available or in processing. 
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(u o H We have thus far been unable to acquire the Informatiou necessary to 

examine the raw geometric properties of either MSS or TM images, nor the raw 
o2 ’^A<^lo4ctric properties of the MSS. Information necessary to examine the low 
2^2 frequency distortions of the Thematic Mapper are not available, e.g., the 
^ M <4 4 B start-, end-, and mid-scan times for the scan mirror are stripped out in 
2 ^ 0 } processing and are not available on A, B, or F tapes. The high 

I frequency spacecraft attitude control information is also not available. For 

oo SS n« m the both the MSS and TM, the internal calibration signals are removed from the 

before the CCTs are made. 


The internal calibration data for the Thematic Mapper is thought to be 
available on CCT-ADDS, for which we have yet to receive any documentation. 
The Goddard group that generates the ADDS tape has told us that this tape is 
considered an "internal" tape, and they are not allowed to send us any 


InfotMtlon *bout It without hlghor lovol Approval • For tha MSS, a faw 
callbratloA raapooaa pointa ara thought to ba in tha trailing part of tha aean 
llna; wa hava not yat auccaaafully Idantlflad thaaa* Alao, wa hava not yat 
baan abla to acqulra a radloaatrlcally uncorractad MSS scana* 

In light of tha data actually availabla froa thaaa two Inatruaanta, our 
propoaal was ovarly aabltloua In tema of axaalnlng tha raw data. Wa will 
thua eoneantrata aora on gaoaatrlcally procaaaad data, using tha uncorractad 
data prlaarlly for radloaatrlc analysat. 

Prograaa aada to data Is auaaarltad in tha abstracts of prasantatlons to 
ba glvan In tlw February aaatlngs (Attachaaant A and B). Additional work, 
underway, but not discussed In tha abstracts Includes tha following: 

Initial aasassaant of band'*to-band raglstaratlon has baan accoapllshad by 
use of color coaposltes and snail area auto-corralatlon tachnlquas. 

Kxaalnatlon of tha spectral equivalency of tha first four bands of tha 
Thematic Mapper with tha four bands of tha Nultlspactral Scanner haa 
baan started. 

Geometric analyses of tha Washington, D.C. scene has begun. A 1:100,000 
scale enlargement of an east-treat strip across the middle of tha scene 
has bean made, and identification of control points begun. The 7 T/2 
minute U.S.G.S. topographic quad for this area have been ordered. 

A generalised routine for examining the contents of the label files and 
non-video data files has been Implemented. Several discrepancies from 
the documentation have been noted and are described In Attachment C. 

A nearest neighbor routine to remove the line length variations In MSS 
Images rapidly, as well as account for all Image distortion accepting 
spacecraft attitude has been developed. 

A scene for which TN and MSS data were acquired simultaneously as well as 
Individually has been identified, and the data ordered. 

Night scenes and day-time ocean scenes required for radiometric purposes 
have been Identified; we plan to use the Landsat assessment system 
(LAL">) to extract the data desired from these scenes. 

The possibility of calibrating the Thematic Mapper using objects other 
than the Earth, and hence avoiding the uncertalnltles of the terrestrial 
atmosphere. Is being examined. 

A formal letter is in preparation modifying our data requirements In 
light of knowledge of the types of data actually available, scenes already 
available, and the availability of scenes In the western U.S. through the 
Canadian receiving station. 

Recomendatlons. Distribute a project status report to each principal 
Investigator on a monthly or more frequent basis; Include a list of data 
availability. Developed the capability to produce rudimentary statistics for 
a scene without requiring full processing; this would allow assessment of the 
utility of a scene for a particular purpose with a much smaller commitment of 
resources. Develope the capability to readily extract and distribute a small 
number of lines of raw data from Thematic Mapper scenes for use In calibration 
and analysis studies. 
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Thsastic Mippsr Intrsbsnd Ksdloastrlc Psrforasncs 
Hu^ H. Klsffsr, Eric Bllsson, Pst Chsvsi, 

Rsyaond Bstson & Wsrrsn Eorgsson 
U. S. Geologic Survey, Plsgstsff, Arisons 86001 

In this ssriy report, «w have concentrated on a rsdloastrlcslly 
corrected, geoMtrlcslIy raw (A data) and a fully processed (P data) laages of 
the Usshlngton D.C* area (ID 4C109-1SU0). In Ardata, every 16th laage line 
la froa a single detector. To allow recognition of features In A data, a 
first order geoaetrlc correction consisting of a A6 saaple shift between 
forward and reverse scan directions was applied. The radloMtrlc correction 
Incorporated in Ardata consists of applying a Radloaetlc look'*up table (&LUT) 
constructed to equalise throughout a scene the aean and standard deviation of 
each detector, with the constraint of nevsr cosipresslng two ON levels Into 
one. The sane RLVT Is used for the forward and reverse scan directions. 

For our statistical study, we tuied a set of lines (512 A, 1024 P) across 
the center of the scene. A sub-scene centered over Chesapeake Bay and 
approxlaatlng a flat field was used for characterising subtle radionetrlc 
differences and noise probleas. 

There are saall differences between the average DN for the 16 detectors 
in each band, typically 0.8 DN; the standard deviations also differ typically 
by 0.4 DN. Differences of node in a band are as large as 6 ON. Soae DN 
levels appear to be strongly favored over adjacent levels. Differences 
between forward and reverse scans are spproxlsMtely 0.1 DN in the wean, and 0 
to 0.6 DN In the standard deviation. Geoaetrlcally resawpled laages (P data). 
In which the fixed correspondence between lines and single detectors Is lost, 
are statistically slalllar to the enseable of detectors In each band In A 
data. 

The Mgnltude of Inter-detector variation Is readily seen by asking an 
laage of the first derivative In the vertical (line) direction of a flat field 
and stretching progressively wider ranges to gray; aost detectors differ froa 
their neighbors by I to 2 DN. 

The effective resolution In radiance Is degraded by * tendency for the TM 
to avoid certain DN levels by about a factor of two. These levels are 
consistent over all bands and detectors, and are spaced by an average of 4 
DN. In band 6, level 127 Is avoided by a factor of 30. This behavior Is 
Basked by resaapllng In the P data. 

At high contrast boundaries, soae of the detectors In Band 5 coaaonly 
over- or under-shoot by several DN and require on the order of SO saaples to 
recover; this behaviours occurs erratically. 

A coherent sinusoidal noise pattern is evident In Detector 1 of Fand 3. 
Qne-dlaenslonal Fourier transforms show that this “stlching'* pattern has a 
period of 13.8 samples with a peak-to-peak amplitude ranging froa 2 to 5 DN. 
Oscillations of the saae frequency, but about half this aaplltude, occur for 
one other detector In this Band. Noise with a period of 3.24 samples Is 
pronounced fwr aost detectors In Band 1 and for 2 detectors In band 3, 
Insignificant In Band 2, weak In Band 4, and obscured by noise In Bands 5 and 
7. 

A set of adjacent detectors In Band 1 has a change of response with a 
period of several scans. These detectors will have a reponse slalllar to 
their neighbors for several scans, then for several scans will have higher or 
lo%ier DN nuabers than their neighbors. There Is no apparent pattern to this 
drift. 




Ar3S Status 1taport» 1933 Fab 18 
Attachaant A 


The randoa noise level of each detector was characterised by the standard 
deviation of the first derivative In the saeple direction across a flat field* 
By this neasure* the noise level 1st Band 1*P to 2*8 DN; Band 2, 0*74 ON 
representative, one detector each at 1*1 and 1*7; Band 3, 0*8 representative, 

4 detectors at 1*0 to 1*3; Band 4, 1*0 representative. Band 5, 1*5 to 2*1; 

Band 7, 1.8 representative, one detector at 3*1* 

The resaepllng algorlthe used In generation of P data Incorporates the 0 
DN values which occupy non-data areas along the left and right edges of a 
fraae* This causes * fussy* fraee edges and adds anoealous low and high DN 
values In the resulting frees histogram* In order to ovoid adverae effects In 
applictlons Involving scene statistics, such as clustering techniques, these 
valutas should be reset to sero after the resampling* 

A principal component analysis Indicates that a composite of the first 
component of Bands 1, 2, and 3; of Band 5 and 7; and Band 4 contain 97Z of the 
Information In the reflectance bands and has reduced the effect of noise* 

A periodic noise removal algorithm developed by Chaves, et* al* has been 
adapted and applied to the coherent noise of Band 3, detector 1; the technique 
appears to work quite well* 

The geometric fidelity of the GSFC fllmirrlter used for Thematic Hspper 
Images was assesed by measurement with accuracy better than 3 micrometers of a 
test grid* The film output has scale errors In both the sample and line 
directions corresponding to 4 pixels In a TM Image, and skew corresponding to 
3 pixels. 
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Landsat-i Multlspactral Scannar Radltmatrle Intraband Parforaanca 
Brlc Bllaaon, Hugh Klaffar 6 Kay Edwards 
(!• S. Gaologlc Survey, Flagstaff, Arisons 86001 

A geoaatrlcally raw laage of Washington D*C. was squired froa BROS Data 
Center. These data have been radloaetrlcally corrected by CSFC; radio- 
aetrlcally ruw data are not available. The data squired show little of the 
detector striping coaaon In earlier )iiltlspectral banner laages. 

The radloaetrlcally corrected data have uniform aeans and standard 
deviations for the detectors in each band; however, the data for different 
detectors utilise a different pattern of DN levels, (typically I DN level out 
of 3 Is not used), resulting in ubiquitous striping of 1 ON amplitude. 

All bands have pronounced coherent noise obvious In contrast enhanced 
Imagea, and evident In natural color composites. One-dlmenslonal Fourier 
transforms of a flat field show a sharp peak with a period of 3.5 samples. 
There Is evidence of weaker period noise at other frequencies. We have 
attempted to remove the coherent noise by both transform techniques and by a 
specially designed sequence of ’’box" filters. 
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DoeuMntation •rrors. 

CBS 10033 ravltlon B 1982 Jun* 11 

k di8cr«p«ncy exists between the docuoMntetlon of the laage Hsjor FrsM Fonet 

(teble 3.5-13) end what actually exists on the digital data* In general the data foraat 

seeas to aatch the docuaentation except as noted below. 

1) Coaputed Cal Leap Value bytes 6281~6284 

This field contains no aeaningful nuaber* Bytes are 
filled with ones giving a floating point nuaber oi 
.33908 B-19 

2) Coaputed Shutter Value bytes 6285-6288 

This field contains no aeaningful nuaber. Bytes are 
filled with ones giving as floating point nuaber of 
.33908 B-19 

3) The fields "Cal Laap Gain Value", "Cal Laap Bias Value", "Applied Gain Value" and 
"Applied Bias Value" do not appear to be located in the byte positions as specified in 
the Table 3*5-13 docuaentation* The probable byte positions are shown below: 

Cal Laap Gain Value - bytes 6293-6296 - docuaentation says byte positions 6289-6292 

Cal Laap Bias Value - bytes 6297-6300 - docuaentation says byte positions 6293-6296 

Applied Gain Value - bytes 6301-6304 - docuaentation says byte positions 6297-6300 

Applied Bias Value - bytes 6305-6308 - docuaentation says byte positions 6301-6304 

4) The byte positions 6289-6292 contain an undocuaented floating point nuaber of unknotm 
significance* This floating point value is roughly 1.00E03 and changes a few percent 
froa one line to the next* 

LANDSAT-D ASSESSMENT SYSTEM LIBRARY COMPUTER COMPATIBLE TAPE (LASLIB-CCT/LAS-CCT) TAPE 
FORMAT DOCUMENTFOR THEMATIC MAPPER. WBS No. 10T031, July 1982 

1) Correction on Table #22 293-298 IMGMAJ 

298 should be 296 

2) On Table #22 Bytes 165-182 SCSTID, SCCNID, SCENID contain no 

aeaningful nuaber when printed with Ll or Al 
foraats* 

3) On Table #22 Bytes 1205-6004 DETDATA, Poraat not specified 

clearly* Values look correct when written out as 
(2R4, 1612). 

4) On Table #17 & #18 Bytes 5-6276 TMHSKl. 

Values look correct when written with a foraat of 
(14, 21R4, 721A1) only they do not aatch any table 
In Ref 2 (GES 10033 Rev B). 



